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Abstract: High-energy neutrinos, arising from decays of mesons produced through the coIUsions of cosmic ray particles 
with air nuclei, form the background in the astrophysical neutrino detection problem. An ambiguity in high-energy 
behavior of pion and especially kaon production cross sections for nucleon-nucleus collisions may affect essentially the 
calculated neutrino flux. We present results of the calculation of the energy spectrum and zenith-angle distribution of 
the muon and electron atmospheric neutrinos in the energy range 10 GeV to 10 PeV. The calculation was performed 
with usage of known hadronic models (QGSJET-II, SIBYLL 2.1, Kimel & Mokhov) for two of the primary spectrum 
parametrizations, by Gaisser & Honda and by Zatsepin & Sokolskaya. The comparison of zenith angle-averaged muon 
neutrino spectrum with the measurement data in IceCube experiment make it clear that even at energies above 100 TeV 
the prompt neutrino contribution is not so apparent because of tangled uncertainties of the strange (kaons) and charm 
(D-mesons) particle production cross sections. An analytic description of calculated neutrino fluxes is presented. 
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1 Introduction 

Atmospheric neutrinos (AN) appear in decays of mesons 
(charged pions, kaons etc.) produced through colUsions of 
high-energy cosmic rays with air nuclei. The AN flux in the 
wide energy range remains flie issue of the interest since 
the low energy AN flux is a research matter in the neu- 
trino osciflations studies, and the high energy atmospheric 
neutrino flux is now appearing as the background noise for 
astrophysical neutrino experiments [1, 2, 3, 4, 5, 6, 7]. 

In spite of numerous AN flux calculations are made (for ex- 
ample [8, 9, 10, 11, 12, 13, 14, 15], see also [16, 17] for a 
review of ID and 3D calculations of the AN flux) there are 
questions concerning to the flux uncertainties originated 
from hadronic interaction models as well as from uncer- 
tainties in primary cosmic ray spectra and composition in 
the "knee" region. 

In this work we present results of the atmospheric neutrino 
flux calculation in the ranee 10-10"^ GeV made with use of 
the hadronic models QGSJET-II 03 [18], SIBYLL 2.1 [19] 
as well as the model by Kimel & Mokhov (KM) [20] that 
were tested also in recent atmospheric muon flux calcu- 
lations [21, 22]. We compute here the differential en- 
ergy spectrum of the conventional neutrinos averaged over 
zenith angles to compare with the data of the Frejus [1], 
AMANDA-II [4] and IceCube [5] experiments. 



2 Calculations vs. the experiment 

The calculation is performed on the basis of the 
method [23] of solution of the hadronic cascade equations 
in the atmosphere, which takes into account non-scaling 
behavior of inclusive particle production cross-sections, the 
rise of total inelastic hadron-nuclei cross-sections, and the 
non-power law primary spectrum (see also [15, 21, 22]). 
Along with major sources of the muon neutrinos, tt^2 and 
K^2 decays, we consider three-particle semileptonic de- 
cays, K^^, K'^^, the contribution originated from decay 
chains K tt ^ {K^g — > tt+tt^, — > tt^tt^), 
as well as small fraction from the muon decays. One can 
neglect the 3D effects in calculations of the atmospheric 
muon neutrino flux near vertical at energies i? > 1 GeV 
and at > 5 GeV in case of directions close to horizontal 
(see [13, 14]). As the primary cosmic ray spectra and com- 
position in wide energy range we use the model recently 
proposed by Zatsepin & Sokolskaya (ZS) [24], which fits 
well the ATIC-2 experiment data [25] and supposedly to be 
valid up to 100 PeV. The ZS proton spectrum at S > 10*^ 
GeV is compatible with KASCADE data [26] as well the 
helium one within the range of the KASCADE spectrum 
obtained with the usage of QGSJET 01 and SIBYLL mod- 
els. Alternatively in the energy range 10 — 10^ GeV we use 
the parameterization by Gaisser, Honda, Lipari and Stanev 
(GH) [17], the version with the high fit to the helium data. 
Note this version is consistent with the data of the KAS- 
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Figure 1: Spectra of the conventional muon and electron 

neutrinos calculated for vertical and horizontal directions. Figure 3: Zenith-angle enhancement of the (j/^ + v^) flux. 



CADE experiment at i?o > 10^ GeV that was obtained 
(through the FAS simulations) with the SIBYLL 2.1. 

Figures 1, 2 display the scale of difference between the con- 
ventional {v^ + v^) spectra and (ve + Ve) one, calculated 
with usage of QGSJET-II, SIBYLL and KM model for GH 
and ZS primary spectra. The difference of neutrino flux 
predictions related to choice of hadronic models is clearly 
apparent. 

Zenith-angle distributions of the conventional neutrinos, 
(j)^^ {E, e) {E, 0°), for the energy range I-IO'^ TeV are 
shown in Fig. 3. Calculations are made with QGSJET-II 
and SiB YLL 2. 1 models both for GH and ZS (ATIC-2) pri- 
mary spectra and composition. As was expected, a shape 
of the angle distribution visibly depends on the neutrino 
energy (at E < 100 TeV) especially for directions close to 
horizontal. The effect of the hadronic models (as well as 
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Figure 2: Conventional neutrino flavor ratios calculated 
with SIBYLL and QGSJET-II hadronic models for the GH 
primary spectrum. 



Figure 4: Conventional (i/^ + D^j) spectrum averaged over 
zenith angles. Curves: the calculation with usage of 
QGSJET-II. Symbols: data of experiments, Frejus [1], 
AMANDA-II [4], IceCube [5]. 

of the primary spectrum) on the angle distribution is rather 
weak, while the spectra differences amount to 80% [15]. 

The calculation of the conventional {v^ + D^) flux av- 
eraged over zenith angles as compared with Frejus [1] 
(squares), AMANDA-II [4] (circles), and IceCube [5] (tri- 
angles) measurement data is shown in Figs. 4, 5. Figure 4 
displays the conventional (i/^ + Dfj) spectrum (averaged 
over zenith angles in the range 0° — 84°) calculated with us- 
age of QGSJET-II model for GH primary spectra and com- 
position (solid line) as well as for ZS one (dashed). The 
difference in neutrino flux predictions resulted from the 
primary cosmic ray spectra becomes apparent at high neu- 
trino energies: the flux obtained with QGSJET-II for GH 
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spectrum at 600 TeV is nearly twice as large as that for ZS 
spectrum. At 1 PeV this discrepancy increases to the factor 
about five. Comparison of QGSJET-II and SIBYLL pre- 
sented in Fig. 5 shows that the former seems more prefer- 
able to desribe the IceCube measurements at the energies 
below 40 TeV (conventional neutrinos). 

The usage of QGSJET-II and SIBYLL models leads to ap- 
parent difference in the neutrino flux, as well as in the case 
of SIBYLL as compared to KM (unlike the muon flux, 
where SIBYLL and KM lead to very similar results [21]). 
On the contrary, the QGSJET-II neutrino flux is very close 
to the KM one: up to 100 TeV the difference does not ex- 
ceed 5% for the GH spectrum and 10% for the ZS one at 
9 = 0°. Note that the muon flux discrepancy in QGSJET-II 
and KM predictions is about 30% at vertical [21]. 

Figure 6 shows the sum of the conventional flux (calculated 
for GH spectrum with usage of QGSJET-II) and prompt 
muon neutrino flux predictions [27] (see also [11, 16, 28]) 
due to nonperturbative models, the recombination quark- 
parton model (RQPM, dotted line) and the quark-gluon 
string model (QGSM, dashed line). The case of ZS spec- 
trum one can see in Fig. 7. The prompt neutrino fluxes were 
obtained [27] with NSU primary spectrum [29], therefore 
they can serve here as upper limits for the prompt neutrino 
flux due to RQPM and QGSM. It worth noting that evalu- 
ation of the prompt neutrino flux obtained with the dipole 
model [30] is close to the QGSM prediction [27] above 1 
PeV. The prompt neutrino flux due to QGSM in the energy 
range 5 TeV < Ei, < 5 ■ 10"^ TeV was approximated by the 
expression 



\2. 011-0.165 



(1) 



where A 1.19 • lO^^** (GeV cm^ s sr)-\ = 100 TeV. 
In this range we neglect the weak angle dependence of the 
prompt neutrino flux. 



10"' 



1 10-= 



10-^ 
10"^ 
10' 
10-^° 




QGSJET+ZS 

QGSJET+ZS+OGSM 

QGSJET+ZS+RQPM 



2 2.5 3 3.5 4 4.5 5 5.5 6 

lg(E/GeV) 



Figure 7: Flux of the conventional and prompt muon neu- 
trinos (case of QGSJET-II + ZS spectrum). 

3 Approximation formula 

Numerical results of the conventional muon neutrino spec- 
tra in the energy range lO'^ — 10^ GeV for different zenith 
angles can be approximated with accuracy (3 8)% by the 
formula: 



5]5]afe„a;V- (2) 



fc=0 ri=0 



Here ^i,{E,y,9) is the flux with units of GeV"^ s"^ sr"^ 
cm^^, X = cos 6*, y = logj^Q(£'ix/GcV). Coefficients akn 
are given in Tables 1,2. 
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Table 1: Coefficients akn for the QGSJET+GH choice 





n 


k 





1 


2 


3 













0.72550 


-6.45625 


4.10284 


-0.87026 


1 


-3.21166 


6.38522 


-3.31293 


0.38300 


2 


1.00337 


-2.46611 


0.99745 


0.01675 


3 


-0.19397 


0.35758 


-0.07515 


-0.04540 


4 


0.01211 


-0.01753 


-0.00135 


0.00537 



Table 2: Coefficients akn for the QGSJETh-ZS choice 





n 


k 





1 


2 


3 













-3.21881 


-7.00088 


4.64475 


-1.07882 


1 


1.60632 


6.92858 


-3.87209 


0.60250 


2 


-1.11835 


-2.65863 


1.20056 


-0.06412 


3 


0.20848 


0.38659 


-0.10641 


-0.03279 


4 


-0.01577 


-0.01909 


0.00037 


0.00467 



4 Summary 

The calculations of the high-energy atmospheric muon 
neutrino flux demonstrate rather weak dependence on the 
primary specrtum models in the energy range 10—10^ GeV. 
However the picture appears less steady because of sizable 
flux differences originated from the models of high-energy 
hadronic interactions. As it can be seen by the example of 
file models QGSJET-II and SIBYLL 2.1, flie major factor 
of the discrepancy in the conventional neutrino flux is the 
kaon production in nucleon-nucleus collisions. 

Comparison of calculated muon neutrino flux with the 
spectrum measured by IceCube shows that QGSJET-II 
is preferable model irrespective of the primary spectrum 
choice. The prompt neutrino contribution due to quark 
qluon string model (QGSM) added to the conventional one 
lead to better agreement with the IceCube measurement 
above 100 TeV. 
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